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Abstract 

Polarimetry with permille-level precision is essential for future electron-positron 
linear colliders. Compton polarimeters can reach negligible statistical uncertainties 
within seconds of measurement time. The dominating systematic uncertainties orig¬ 
inate from the response and alignment of the detector which records the Compton 
scattered electrons. The robust baseline technology for the Compton polarimeters 
foreseen at future linear colliders is based on an array of gas Cherenkov detectors 
read out by photomultipliers. In this paper, we will present a calibration method 
which promises to monitor nonlinearities in the response of such a detector at the 
level of a few permille. This method has been implemented in an LED-based cali¬ 
bration system which matches the existing prototype detector. The performance of 
this calibration system is sufficient to control the corresponding contribution to the 
total uncertainty on the extracted polarisation to better than 0.1%. 


1 Introduction 


Since the event rates of electroweak processes depend on the chirality of the colliding par¬ 
ticles, beam polarisation is a key ingredient of the physics programme of future electron- 
positron colliders [1] and its precise knowledge is as important as the knowledge of the 
luminosity. Thereby, the luminosity-weighted average polarisations at the interaction 
point {V)inrau which is the value relevant for the interpretation of collider data, have to 
be determined by combining fast measurements of Compton polarimeters with long-term 
scale calibration obtained from reference processes in collision data. 

In particular for the International Linear Collider [2], where both beams are foreseen to 
be longitudinally polarised, it is required to control (P)iumi with per mille-level precision 
in order not to limit the precision of cross-section measurements. Two Compton po¬ 
larimeters [3] per beam aim to measure the longitudinal polarisation before and after the 
collision with a precision of 8VjV < O.259{0. It should be noted, though, that this goal 
is not driven by physics requirements, but by what used to be considered feasible experi¬ 
mentally. Thus, further improvements in polarimetry would still have direct benefits for 
the physics potential of the machine. 

In order to evaluate the effects on the mean polarisation vector caused by the beamline 
magnets between the polarimeters and the interaction point, by the detector magnets and 
by the beam-beam interaction, spin tracking simulations are required. These effects have 
recently been studied in [1], concluding that not only a cross-calibration of the polarime¬ 
ters to 0.1 % is feasible, but also individual extrapolations of upstream and downstream 
measurements to the e’''e“ interaction point. 

From reference reactions, the long-term average of the polarisation at the interaction 
point can be determined. In particular W pair production and single W produc¬ 
tion [S] have been studied in this respect, showing that precisions of about 0.15% can be 
achieved after several years of data taking. These results are quite robust with regard 
to systematic uncertainties: e.g. in the measurement using W pair production, even for 
a conservative assumption of 0.5% uncertainty on the selection efficiency for the signal 
and 5 % for the background, the impact on the uncertainty of the polarisation measure¬ 
ment was found to be below the statistical uncertainty for an integrated luminosity of 
500 fb“^ [6j. Any imperfection in the beam helicity reversal, i.e. differences in the mag¬ 
nitude of the polarisation between measurements with left- and right-handed polarised 
beams, has to be corrected for based on the polarimeter measurements. 

The two polarimeters per beam are located about 1.8 km upstream and 160 m downstream 
of the interaction point, providing non-destructive measurements of the longitudinal beam 
polarisation based on the polarisation dependence of Compton scattering. They have been 
designed for operation at beam energies between 45 GeV and 500 GeV. Several options for 
detecting the Gompton-scattered electron^ are being considered. The baseline solution is 
a gas Gherenkov detector consisting of 20 channels, each 1 cm wide. A two-channel proto¬ 
type of such a detector has been successfully operated in testbeam, where the alignment 

^for typical ILC beam polarisation values of Vf.- > 80% and 'PeS- > 30% or even > 60%. 

^or positrons in case of the positron beam of the ILC which is equipped analogously. Throughout this 
paper, we’ll use “electron” to refer to both cases, unless explicitly stated otherwise. 
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requirement for the ILC has nearly been reached and electronic noise was found to be at 
a negligible level [9]. The remaining challenge for per mille-level polarimetry with a gas 
Cherenkov detector is the linearity of the device. 

In this paper, we present a calibration system which is able to monitor nonlinearities with 
sufficient precision and provides the abilities to correct for them. The paper is organised 
in three main parts: In section [2l we summarise the most important facts about the ILC 
Compton polarimeters before introducing the baseline detector for the Compton scattered 
electrons and discussing the systematic uncertainties of the measurement, which define 
the target precision of the calibration system. In section El we introduce a suitable 
linearisation method, discuss the resulting requirements, and present a realisation of a 
suitable light source. The performance of the calibration system and its application to 
the ILC Compton polarimeters will be discussed in sectionlUbefore concluding in section|5l 


2 Compton Polarimetry at the ILC 


Compton scattering of circularly polarised laser light off the lepton beam is the method 
of choice for polarimetry at high-energy colliders for a number of reasons: As a pure QED 
process, the scattering cross-section can be calculated precisely, with radiative corrections 
less than 0.1% m, and offers large polarisation asymmetries. The interaction with the 
laser conserves the beam quality so that polarimetry does not interfere with collision data 
taking. The laser polarisation can be flipped very quickly using e.g. a Pockel’s cell and its 
absolute value can be controlled at the level of 0.1% m- Backgrounds will be determined 
in-situ by shutting the laser off and thus can easily be corrected for. 

Compton polarimetry was successfully employed at many previous accelerators, for in¬ 
stance at the SLC [H] or at HERA [T2l|T3]. In this section, we will summarise the most 
important generic aspects of Compton polarimetry following the notation from [H], be¬ 
fore summarising the current design of the polarimeters foreseen at the ILC. Additional 
details on the ILC polarimeters can be found in |3]. 

2.1 Compton Polarimetry Basics 

The kinematics of Compton scattering are fully determined by two variables, e.g. the 
azimuthal and polar angle of the scattered electron. In particular, the polar angle is 
equivalent to the energy of the scattered electron. Besides the unpolarised part, the 
double-differential Compton cross section na contains terms proportional to the linear 
light polarisation and to the product of circular light polarisation and transverse or longi¬ 
tudinal electron polarisation. When integrating the cross section over the full azimuthal 
angle, only the term containing the longitudinal polarisation remains, so that any sensi¬ 
tivity to linear light polarisation and transverse electron polarisation is removed. Thus 
transverse beam polarisation requires a double differential measurement in energy and 
e.g. vertical positionH, while for longitudinal polarimetry, which is the primary concern 

^For a recent discussion of transverse polarimetry at the ILC c. f. m- 
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at the ILC, it is sufficient to measure either the energy of the scattered electron or the 
scattered photon. 

The minimum energy of the electrons is given by the Compton edge energy 

-E'min = - ) ( 1 ) 

1 + X 


where Eq is the beam energy, while the dimensionless variable x relates Eq to the energy of 
the laser photons wq, the electron mass m and the crossing angle 9q between the electron 
beam and the laser: 


X = 






cos^ {9q/2) ~ 






( 2 ) 


For e.g. a green laser (cuo = 2.33 eV) shining head-on to an electron beam with Eq = 
250 GeV, X ~ 9 and thus i^min — O.li^o = 25 GeV, whereas the remaining energy is 
carried away by the photon. 


In the laboratory frame, both the scattered electrons and photons are contained in a 
narrow forward cone around the electron beam direction, with typical opening angles of 
10 — 20/rrad. 

The differential Gompton scattering cross section with respect to the normalised photon 
energy y = oo/Eq can be written as [T5] 


da{X,Vz) 2(To 


+ 1 — y — 4r(l — r) + VzXrx{l — 2r)(2 — y) 


dy X \_1 — y 

or with respect to the absolute energy of the Gompton scattered electron it becomes 

dcr 1 dcr 
dE Eq dy ’ 

where 

y 


(Jo = 249.5 mb and r = 


x(l - y)' 


( 3 ) 


( 4 ) 


( 5 ) 


Vz is the longitudinal electron beam polarisation and A describes the laser polarisa¬ 
tion. Thus, the differential cross section is separated into a spin-independent and spin- 
dependent part. For convenience of the reader we note that r = 1 corresponds to the 
Gompton edge. 

From equation [3] it can be seen that the asymmetry between the differential cross sections 
corresponding to the two laser helicity states sgnA = ±1 is proportional to the product 
of laser and beam polarisation: 


A{Vz) 


da-(PQ-da+{VQ 
da-(V,) + da+{V,) ' ' 


( 6 ) 


Herein, ((j“) denotes the case of parallel (anti-parallel) laser and electron helicities. 
The proportionality factor is given by the inverse of the analysing power AV, which is 
the asymmetry expected for XVz = 100 %. 
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Figure 1: (nj) Differential Compton scattering cross section for a beam energy of Eq = 
250 GeV and green laser light. The different colonrs show the dependence of 
the differential cross section on the prodnct of laser and beam polarisation. 
dbl) Asymmetry for different absolnte values of the beam polarisation. The 
asymmetry corresponding to \Vz ■ X\ = 1.0 (black) is also called analysing 
power. 


Figured] shows the differential Compton scattering cross section as well as the correspond¬ 
ing asymmetry for different beam polarisations with respect to the energy of the Compton 
scattered electron. In these hgures, the beam energy is set to Eq = 250 GeV and a green 
laser is chosen, as this reflects the ILC baseline design parameters. 

In practice, da will be determined by measuring the count rate in a detector channel 
with hnite size, typically after converting the energy spectrum of the Compton scattered 
electrons into a position distribution via a magnetic chicane, as will be explained in the 
next section. Per channel i, the observed asymmetry is given by 




N-(V.) - N*iVQ 
N-(V,) + N*(V,) 


(7) 


assuming the same luminosities for both laser helicities. 


The measured polarisation values from the individual channels can then be combined 
into a weighted average, where the weights Wi are chosen to minimise the total statistical 
uncertainty, which is achieved by giving the largest weights to the channels with largest 
analysing power AVi m- 


1 

{Vz-AVi)-^-l 


( 8 ) 
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2.2 The ILC Compton Polarimeters 


The ILC Compton polarimeters will be operated in a multi-event mode with 0(10^) 
Compton interactions per bunch crossing. For the polarisation measurement, both the 
scattered photon and electron can be analysed equivalently. At HERA, the longitudinal 
polarimeter operated successfully in a multi-photon mode by detecting the sum of the 
photon energies in a calorimeter na. At ILC energies and intensities, however, the 
cumulative energy of 0(10^) photons carrying on average (9(10) GeV each would amount 
to 0(10) TeV per bunch crossing, all within a narrow cone of 0(10) /rrad. Therefore, it 
seems highly unlikely that photon detection could contribute to precision polarimetry at 
the ILC. 

The situation is very different for the detection of the scattered electrons as employed at 
SLD. As charged particles, the electrons can be energy-analysed by dipole magnets. This 
translates the energy measurement into a position measurement which does not require 
a total absorption of the energy. Thus, the shape of the spectrum can be monitored 
under multi-event conditions. Additionally, regions of high or low analysing power can 
be evaluated separately and enter in the total measurement with an appropriate weight 
(c.f. equation [8]). 



Figure 2: Sketch of the magnetic spectrometer for the ILC upstream polarimeter 
(from [5]). The Compton-scattered electrons with the lowest energies are 
deflected most by the spectrometer. 

The polarimeter chicanes consist of a sets of dipole magnets. Figure |2] shows as an example 
a sketch of the chicane foreseen for the ILC’s upstream polarimeters. The magnetic held 
strengths of the magnets are chosen such that the fan of Compton scattered electrons are 
dehected far enough from the main beam so that a segmented counting detector can be 
safely placed and operated next to the beam pipe. Both the upstream and the downstream 
polarimeter chicanes will be operated at a constant magnetic held for all beam energies. 
This implies that the Compton interaction point moves laterally with the beam energy, 
but ohers the advantage that the position and size of the Compton spectrum in the 
detection plane stays invariant. Thus, the same measurement quality can be maintained 
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for all beam energies. A detailed description of the polarimeter chicane designs for the 
up- and downstream polarimeter is discussed in reference [3]. 


2.2.1 The Upstream Polarimeters at the ILC 

At the upstream polarimeter, the beam spot of the yet undisturbed bunches is small and 
low background levels are expected [H] . Therefore, conventional lasers featuring moderate 
pulse energies can be used in order to probe every single bunch in an ILC bunch train. 
Different ILC beam parameter sets foresee between 1312 and 2625 bunches per bunch 
train, where the trains last 1ms each and occur with a repetition rate of 5 Hz [2]. The 
ability to probe every bunch in a bunch train allows to accumulate high statistics in a very 
short time. Furthermore, it allows to resolve possible variations in the bunch polarisation 
within a bunch train which could be present due to beam dynamics. 

On the downside, the expected radiation exposure of the detector from the Compton elec¬ 
trons alone amounts to up to about 1 MGy per year, which requires to employ sufficiently 
radiation hard detector technologies. The resulting readout rate in the order of 10 MHz 
requires a fast responding detector as well as fast readout electronics. 

The upstream polarimeter chicane as sketched in hgure [2] spreads the Compton spectrum 
over approx. 20 cm in the horizontal plane. Due to the benign conditions, the instrumen¬ 
tation can cover nearly the whole spectrum down to about 2 cm from the undisturbed 
main beam, corresponding to Compton electrons up to U = 125 GeV. This clearly in¬ 
cludes the zero crossing point of the asymmetry (c.f. figure [lb]), which together with the 
Gompton edge position provides important calibration reference. 


2.2.2 The Downstream Polarimeters at the ILC 

The downstream polarimeter chicane starts 120 m behind the interaction point, 

where the Gompton interaction point is located at a secondary focus point of the machine 
optics. However, due to the energy loss during collisions, the beam cross section is four 
orders of magnitudes larger than at the upstream polarimeter [1] . In addition, significant 
beam background at the level of C(IO^) photons and 0(10^) charged particles per bunch 
crossing is expected [18]. For this reason, the laser power needs to be increased by more 
than three orders of magnitude to still maintain a reasonable signal to background ratio. 
This, in turn, reduces the maximal laser repetition rate to 0(1) laser pulse per ILG bunch 
train. Therefore, the detector readout speed is of no concern. 

For the same reason, the radiation dose due to the Gompton electrons themselves is three 
orders of magnitude lower than at the upstream polarimeter. On the other hand, the 
substantial backgrounds, which are present for every bunch independently of the presence 
of a laser pulse, will add to the total radiation dose. Thus it seems prudent to keep similar 
radiation tolerance requirements as for the upstream polarimeter. 

Since the extraction line is designed with an aperture of 0.75 mrad with respect to the 
interaction point, the polarimeter detector needs to be placed in a distance of 15 cm to 
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the main beam trajectory. In order to still be able to probe a significant part of the 
fan of Compton scattered electrons, the polarimeter chicane is slightly modified and two 
additional magnets operating at a higher field strength are inserted in the spectrometer 
part resulting in a larger dispersion [12] . At a beam energy oi E = 250 GeV, Compton 
scattered electrons from E = AA GeV down to the Compton edge at E = 2b GeV can be 
detected in this setup. The accessible spectrum is spread over a distance of approximately 
13 cm. Due to the expected background, detection methods insensitive to photons and 
with an intrinsic energy threshold for charged particles are considered advantageous. 


2.3 Detectors for the ILC Compton Polarimeters 

A natural choice for fast counting detectors are Cherenkov detectors read-out by a photo¬ 
multiplier. For relativistic electrons, this detection mechanism is intrinsically linear, since 
the amount of emitted Cherenkov light is proportional to the number of electrons passing 
through that channel. Depending on the choice of Cherenkov medium, the threshold for 
Cherenkov radiation can be in the MeV regime, which suppresses backgrounds from low 
energetic electron-positron pairs from beamstrahlung. This is in particular true for gases 
with refractive indices only slightly above unity. E.g. for C 4 ,Eiq with n = 1.0014, the 
Cherenkov threshold is at 10 MeV. Due to their reliability, robustness and proven perfor¬ 
mance in the so far most precise Compton polarimeter which measured the polarisation 
of the electron beam at SLC mi, gas Cherenkov detectors are considered the default 
solution for the ILC polarimeters. As will be discussed below, the detector linearity is 
one of the most important systematic uncertainties of this type of detector. Therefore, 
we present in this paper a calibration system which is suitable to detect and correct for 
nonlinearities at the sub-percent level in-situ. 

Recently, a Cherenkov detector for polarimetry based on quartz has been proposed [20] . 
With a much higher light yield (due to the higher refractive index) in combination with a 
higher granularity, such a system would robust against nonlinearities up to the percent- 
level. However, this concept still awaits a full demonstration in testbeam and a full 
evaluation of the robustness against backgrounds, since the Cherenkov threshold is much 
lower. 

Ultimate granularity would be achieved with a silicon pixel detector, eliminating the 
question of linearity completely and increasing significantly the possibilities to control 
the alignment. However, in this case R&D is still required to reach a proof-of-concept 
level, where the main concern is the high local data rate in the most intense region of the 
Compton spectrum, followed by radiation issues. 


2.3.1 A Gas Cherenkov Detector for Compton Polarimetry 

A prototype of a gas Cherenkov detector for ILC polarimetry has been designed, built 
and operated successfully in testbeam [2]. Based on the experience gained at SLC, the 
design of the channels has been improved to U-shaped aluminium tubes with a square 
cross section of 1 x 1 cm. A schematic side view of one channel is depicted in Figure [3al 
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When the electrons traverse the horizontal tube of the Cherenkov counter, they induce 
Cherenkov light. Mirrors at the end of the tube reflect the Cherenkov light upwards such 
that the actual photo detector can be placed outside the Compton scattering plane. The 
other upward pointing “leg” of the U-shaped channels is foreseen to host an LED-based 
calibration system. For the upstream polarimeter, 20 such channels would be staggered 
along the shallow angle exit window. 



Figure 3: (Hj) Schematic drawing of one gas Cherenkov detector channel, (jbj) Picture of 
the two-channel Gas Cherenkov detector prototype. Both from [9]. 

Figure |3b] shows a drawing of the prototype which consists out of two non-staggered 
channels. In testbeam, it was demonstrated that the aimed detector alignment goals can 
be achieved, in particular when multi-anode photomultipliers (MAPMs) are used which 
allow to measure asymmetries in light distribution in each channel. Since the light pattern 
on the detector plane is related to the detector-beam alignment, the alignment can also 
be measured and monitored during a physics run without spending dedicated beam time 
on alignment runs [9]. 

Figure H] shows the number of expected Compton scattered electrons per polarimeter 
channel for a simulated polarisation oiV = 80%, as obtained from the fast Linear Collider 
Polarimeter Simulation LCPolMC HIEI]. The axis on the right-hand side illustrates the 
expected QDC signal assuming the parameters of the prototype detector [9], thus a light 
yield oir] = 6.5 photoelectrons/ Compton electron, a photomultiplier gain oi g = 3.0 x 10^ 
and a QDC resolution of 25 fC/QDC count. 


2.3.2 Sources of Systematic Uncertainties 

In the multi-event mode, the statistical uncertainty of the polarisation measurement drops 
rapidly such that already after a few seconds of data taking, the polarisation measurement 
uncertainty is dominated by systematic uncertainties. Table [1] summarises the envisaged 
error budget at the ILC in comparison to the obtained systematic uncertainties at the 
SLC polarimeter [m. 
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Figure 4: Simulated number of Compton scattered electrons per polarimeter channel for 
the different laser helicities. The right axis shows the expected QDC signal 
assuming the parameters of the prototype detector, see text. 
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source of uncertainty 

SLC 

ILC goals 

laser polarisation 

0.1% 

0.1% 

detector alignment 

0.4% 

0.15-0.2% 

detector linearity 

0.2% 

0.1% 

electronic noise and beam jitter 0.2% 

0.05% 

Total 

0.5% 

0.25% 


Table 1: Error budget of the polarisation measurement at the ILC [1] in comparison to 
the determined systematic uncertainties at the SLC polarimeter m- 


Laser polarisation. The laser polarisation and its flip accuracy have been controlled to 
the level of 0.1% already at the SLC polarimeter. Currently, no improvement is 
envisaged for the ILC case. 

Noise &; jitter. Contributions originating from electronic noise, which were sizeable at 
the time of the SLC polarimeter, are expected to be negligible for modern DAQ 
hardware as applied in the ILC polarimeters, which is conhrmed by the testbeam 
experience with the prototype |H]. Also the beam jitter is expected to be signihcantly 
reduced in the ILC case, since the luminosity goals require to limit the jitter to at 
most 10% of the beam sizes [2]. 

Analysing power. At cross-section level, the analysing power can be calculated with 
high precision from the Compton cross section known currently to the order [10] . 
Therefore, the knowledge of the analysing power is dominated by the alignment of 
the detector with respect to the Compton spectrum, which includes lateral shifts 
of the Compton IP as well as the hnite knowledge of the magnetic helds of the 
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chicane. Both effects are included when the alignment is derived from measuring the 
unpolarised Compton spectrum, e.g. from the sum of data taken with the two laser 
helicities. The foreseen error budget translates into a vertical alignment precision of 
at least 100 /rm and a precision of the tilt angle of the detector of 1 mrad with respect 
to the design beam orbit |2U]. A granularity of about 1cm provides sufficient detail 
to achieve this level of control by monitoring the shape of Compton spectrum [9]. 
It has been shown that the requirements for shifts as well as tilts of the detector 
could be met with the prototype, in particular when assisted by information of the 
intra-channel light pattern read-out by multi-anode PMTs [9]. 

Detector nonlinearities. Apart from the alignment, the detector nonlinearities are the 
other source of uncertainty which needs to be improved significantly in order to 
meet the precision goal for the ILC. From equation ([7j), it becomes clear that a 
polarimeter detector does not need to be calibrated on an absolute scale, since any 
constant calibration factor cancels. However, the linearity of the detector response 
is essential for per mille-level polarimetry. 


Impact of nonlinearity on polarisation measurement 



Figure 5: Systematic uncertainty on polarisation measurement due to different detector 
nonlinearities. The error band shows the RMS interval obtained from 200 ran¬ 
domly chosen detector response functions. The exact dehnition of nonlinearity 
employed here is described in section 13.11 

In order to investigate the maximal tolerable detector nonlinearity which still keeps 
the nonlinearity contribution to the overall systematic uncertainty below 0.1%, the 
impact of different nonlinear detector response (or transfer) functions on the mea¬ 
surable polarisation has been studied in Monte Carlo simulations. For each value 
of the relative nonlinearity, 200 third order polynomials with random parameters 
have been generated and applied as transfer functions in the simulated polarisation 
measurement. The results are depicted in Figure El The red line indicates the ab¬ 
solute value of the mean deviation between the reconstructed and true polarisation. 
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while the yellow band shows the RMS of those samples. It becomes obvious that 
the detector nonlinearity must be controlled at the level of 0.5% in order to keep the 
resulting contribution to 5VjV below 0.1 %. A method to achieve this is discussed 
in the next chapter. 


3 A Calibration System for the ILC Polarimeters 


In order to control the detector nonlinearity on a level of below 0.5%, novel calibration 
techniques are needed. In this section, we present a calibration system, which is capable 
to meet these requirements. 


3.1 Definitions of Nonlinearity 

Basically, there are two ways of describing the nonlinearity of a detector, which are illus¬ 
trated in Figure El The integrated nonlinearity (INL) specifies the absolute deviation of 
the detector response T(a:) from a linear response L(a:), where x is an applied signal (c.f. 
Figure l6a| . Thereby, the linear response is dehned by the response to a signal Xref which 
marks the end of the range to be calibrated. 

Another way of describing a nonlinearity is the differential nonlinearity (DNL). Here, the 
deviation between the slope of the detector transfer function and the linear response is 
used (c.f. Figure l6b|) . Although the DNL does not preserve any information about the 
absolute scale of the detector response, all the information regarding the linearity of the 
response is still present. Thus, also a DNL measurement can be used to linearise the 
detector response. 




(a) Integrated nonlinearity (b) Differential nonlinearity 

Figure 6: Different ways of specifying the nonlinearity of a detector response. 
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3.2 Permille-Level Linearisation 


Measuring the nonlinearity of a light detector to permille-level is very challenging as 
this usually requires an intensity calibrated light source of the same precision. Different 
calibration methods based on LED light pulses were studied in Reference izg. In the 
following, we will focus on a double pulse method, which pursues a differential calibration 
approach: One light pulse (base pulse x) is tunable in its light intensity over the expected 
dynamic range of the detector. A second light pulse (differential pulse Ax) features a 
constant intensity, which is small compared to the dynamic range. Measuring the detector 
response to x and x + Ax, allows to estimate the DNL. This approach is advantages as 
it does not require a calibrated light source at all as it will be demonstrated in the next 
section. Details on this method are also discussed in Reference [23] . 


3.2.1 Mathematical Foundations and Linearisation Method 

In the following, we outline the mathematical foundations and the procedure of permille- 
level linearisation using the double pulse method. 

The most general transfer function of a detector can be parametrised as 

T(x) = A-I-(R-I-n/(x)) X. (9) 

Herein, A denotes an offset (e.g. due to dark current), which can be easily measured 
and subtracted. Thus, we set A = 0 for simplicity. nl{x) is a arbitrary function which 
describes the nonlinear part of the detector response. In order to linearise the detector 
response, the correction function 


Corr(x) 


B 

B + nl{x) 


( 10 ) 


needs to be determined. 

The derivative of T(x) is proportional to the difference between the values of the detector 
response for a dedicated light signal x and the signal which is increased by a constant and 
small additional quantity of light from the differential light pulse Ax: 

D(x) := T(x) oc A T(x) = T(x -I- Ax) — T(x) (11) 

LliL L^Jb 

If nl{x) is small, which is intrinsically the case for photomultipliers, it holds 

D(x)« D(T(x))0CAT(T(x)). ( 12 ) 

Thus, measuring AT in dependence of T(x) contains (up to an proportionality factor) 
approximately the same information as the direct derivative of the transfer function in 
dependence of the absolute amount of injected light x. Thus, we achieve an approxi¬ 
mate determination of the derivative of the transfer function which does not rely on the 
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knowledge of the injected quantity of light x, but only on the corresponding detector 
responses. 

The measurement of the data pairs (T(x), A T( T(a:))) can now be used to build the 
correction function Corr(x). For this purpose, the data points are htted with a higher 
order polynomial in order to obtain a continuous function which can then be processed 
further. 

Combining Equations ([9]) and flTT]) gives 

A T(x) Ax • 4— T(x) = Ax • (B + nlix) + nl'(x)x). (13) 

dx 

Thus, the average of A T over the whole dynamic range can be identihed as 

(A T(x)) = Ax • {B + (n/(x) + nl\x)x)) = Ax • B. (14) 

Note that the average over the n/(x)-dependent terms in Equation fflTll vanishes, as by 
dehnition n/(xref) = 0: 


- / inlix) + nl'(x)x) dx =-[n/(x)x]o™'^ = 0. (15) 

Xref Jo Xref 

In a next step. Equation ffldll is corrected for the constant part 

A T(x) — Ax ■ B = A T(x) = Ax(n/(x) + nl'{x)x) (16) 

and AT(x) is integrated in order to obtain an expression for nl{x): 

AT(x')dx'= Ax • nZ(x) (17) 

Finally, Equation flTT)l and CHI can be combined in order to construct the correction 
function of Equation ffTOjl . The remaining proportionality factor Ax in fflTll and ffTTH 
cancels. It should be stressed once more that Corr(x) ~ Corr(T(x)) can be determined 
without knowing the absolute input light signal. The only requirement is that the ad¬ 
ditional quantity of light Ax is constant and small compared to the quantity of light x 
which scans the dynamic range. 



3.2.2 Monte-Carlo Study 

In order to test the presented linearisation method, a Monte-Carlo simulation has been 
developed, which takes into account the simulation of the light emission, attenuation, 
detection, and digitisation. The simulation parameters have been taken from datasheets 
of the foreseen DAQ equipment or, in case of the light source parameters, have been 
adjusted in agreement with measured distributions. With this simulation, the performance 
and limitations of the described linearisation method can systematically be studied and 
parameters for a measurement setup can be optimised. 
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As the desired linearisation is in the sub-percent range, a sufficiently large number of data 
points (T(a;), A T( T(a;))) with a sufficiently large statistic needs to be acquired. From 
the simulation we hnd that the number of scan steps within the dynamic range plays a 
small role as long as it exceeds four scan steps. This can be understood since nonlinear 
contributions of higher orders in the transfer function play a decreasing role and, thus, 
more data points do not strongly constrain the polynomial fit anymore. However, in order 
to ensure a good modelling of the nonlinearity by the fit it is advantages to adjust the 
number of scan steps according to the time budget of a nonlinearity measurement. A 
much stronger influence on the performance of the linearisation method has the number 
of measurements taken at one scan step since this effects the statistical uncertainty of 
the individual data points. It has been found that in order to allow for a sub-percent 
linearisation at least 3 x 10® measurements per scan step are desirable. 

One key requirement is the stability of the differential light pulse. Figure [7] shows the 
impact of a drift of the intensity of the differential pulse by 0.2% and 2.0% within the 
scan range on the resulting nonlinearity before and after the linearisation. This figure 
illustrates that for small drifts the detector response can still be well linearised. However, 
once the drift gets too large the linearisation method does not work satisfactory anymore. 
In order to estimate the maximal tolerable drift of the differential pulse, this parameter 
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Figure 7: Simulation of the effect of 0.2% and 2.0% drift in the light intensity of the 
differential pulse during the scan of the dynamic range. The horizontal line 
indicates the design goal of a nonlinearity smaller than 0.5%. For a too large 
drift, the method does not sufficiently linearise the detector response anymore. 


has been scanned in the simulation. Thereby, for each scan step, 200 randomly generated 
nonlinear transfer functions featuring the same total integrated nonlinearity of 0.8% but 
different shapes have been simulated. Figure [H] shows the average effective nonlinearity, 
defined as INLmax — INL^in in the studied range, after the linearisation. The uncertainty 
band indicates the RMS of the resulting distribution around the central value. From the 
figure it gets clear that the differential pulse must not drift by more than 0.5% during 
the scan of the dynamic range in order to still fulhl the goal of a linearisation better than 
0.5%. 
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Figure 8: Effective nonlinearity in dependence of different simulated drifts of the differ¬ 
ential pulse. The drawn error band shows the RMS around the mean value 
from all 200 simulated shapes of nonlinearities. In order to keep the resulting 
nonlinearity below 0.5%, the differential pulse must not drift by more than 
0.5%. 


The presented linearisation method is based on the assumption of only small detector 
nonlinearities (cf. Equation flT2|) h Therefore, also the limitations of the method with 
respect to the underlying nonlinearity needs to be studied. Figure |9] shows the effective 
nonlinearity before and after linearisation for different simulated nonlinearities. One can 
see that nonlinearities of up to 3.5% can still be corrected successfully with the proposed 
approach such that the effective nonlinearity stays below 0.5%. This range is sufficient 
for photomultipliers, which are intrinsically very linear devices. 


3.3 Requirements for the calibration system 

The differential calibration approach has the advantage that the requirements for the 

calibration light source is reduced. However, it still needs to fulfil a few specihcations: 

General design. The calibration system must consist of two independent light sources, 
which can be independently adjusted in their light intensity. 

Stability of differential pnlse. The light pulse of the differential pulse dehnes the mea¬ 
sure of the derivative and, therefore, it must be very stable. In the previous section 
it has been shown that a relative drift of the differential pulse of 0.5% is tolerable 
in order to still meet the calibration requirements. 

Dynamic range of base signal. The light source providing the base pulses has to be 
adjustable in its light intensity over the whole dynamic range of the planned Comp¬ 
ton polarimeters in order to be able to probe the differential nonlinearity within this 
range. As can be seen from figure IU up to 200 Compton electrons per bunch crossing 
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Different simulated non-linearities 



Figure 9: Effective nonlinearity in dependence of the simulated relative INL. The drawn 
error band shows the RMS around the mean value from all 200 simulated 
shapes of nonlinearities. The proposed linearisation method still works for a 
relative INL of up to 3.5%. 


can be expected in a single detector channel, assuming the nominal instantaneous 
Compton luminosity. 

Additionally, there are some more requirements, which are important in order to meet 

the conditions of the actual polarisation measurement at the ILC. 

Wavelength. The light sensors at the Compton polarimeters are supposed to detect 
Cherenkov light, whose wavelength intensity distribution peaks in the UV range. 
Therefore, the calibration light source should emit light in this wavelength range as 

well. 

Short pulses. An ILC bunch with a longitudinal bunch length of 300 /xm [2] traversing 
the Compton polarimeters creates a Cherenkov light pulse of only t = 1 fs. This is 
much shorter than the usual transit time of the electron avalanche in a photomul¬ 
tiplier. Therefore, the light pulse of the calibration source should be as short as 
possible, but at least as short as the photomultiplier transit time, which is of the 
order of 10 ns for photomultipliers we consider in out setups. 

Applicable in detector design. The light source should ideally be very small such that 
it can be integrated in the existing polarimeter design. 


3.4 LED driver 

According to these specifications an LED driver has been developed based on the calibra¬ 
tion light source of the CALICE tile hadron calorimeter [2l]. Figure fTOal shows a picture 
of the LED driver. It is equipped with two UV-LEDs which feature a peak intensity at 
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A = 395 nm and a rather large spectral width of a few 10 nm. The LED light intensity is 
dehned by two external reference voltages 17 ledi, 2 - Both LEDs receive the same trigger 
signal. 
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Figure 10: ^ Picture of the LED driver. The design is chosen such that it is applicable 
in the gas Cherenkov detector design, (jbj) LED driver circuit. 


It is eye-caching that the dimensions of the board are rather small (50 x 8.5 mm^) such 
that the driver can easily be integrated in the existing Cherenkov detector design (see 
section [2.3. ip . Both UV-LEDs are located on the bottom side of the board, whereas the 
remaining electronic components are placed on top. As LEDs are rather temperature 
dependent devices, two temperature sensors are located on the board. An analogue tem¬ 
perature sensitive resistor (PtlOOO) is placed directly opposite the LEDs on the top side. 
A second, digital temperature sensor (DS18B20U) is present at the other end of the board, 
which provides a 12-bit digitised temperature, read out via a 1-wire bus. The board is 
connected via a 10-pin connector. 

A simplified schematic of the LED driver circuit is shown in Figure ITObl This design [21] 
allows for very short light pulses. There are three phases in the generation of a light 
pulse: Firstly, the capacitor C gets charged via the resistors R and Rp {R 3> Rp) driven 
by the external voltage Dled- Secondly, when the trigger arrives, the transistor is closed 
and the capacitor discharges via Rp and i?LED- The LED is emitting light. Thirdly, the 
resistance of the LED rises as the voltage at the capacitor drops due to the nonlinear 
I — V behaviour of the LED. The discharging continues mainly via the parallel resistor 
Rp and, thus, the LED light pulse is quenched. 

The LED board has been carefully characterised in a dedicated test setup. Figure [11] shows 
a schematic overview of the setup. It consists of a light tight box with an optical bench 
inside, on which the LED driver as well as a photomultiplier is mounted. As photodetector 
a segmented photomultiplier by Hamamatsu (R5900-03-iyi4) has been selected. This type 
has also been operated in the gas Cherenkov detector prototype described in section 12.3.11 
For a coarse adjustment of the light intensity, different neutral density Liters can be 
inserted into the optical path. In the standard setup, the photomultiplier is read out by 
a 12-bit QDC. 
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Figure 11: Schematic picture of test setup for characterising the LED driver board and 
for measuring the differential nonlinearity of a photomultiplier. 

In order to resolve the pulse length of light pulses emitted by the LED driver, the photo¬ 
multiplier has been connected to an oscilloscope at Erst. Figure [T2a] shows the recorded 
wave forms for different reference voltages Dled averaged over 5120 individual light flashes. 
It is visible that the measurable pulse length is below 10 ns and, thus, the driver fulhls 
the described pulse length requirement. The observable wiggles at t > 40 ns can be at¬ 
tributed to photomultiplier after pulses and signal reflections because of imperfect signal 
termination at the oscilloscope. In order to exclude those nonlinear effects from the QDC 
measurements it is important to chose the gate carefully. 

In Figure [T2bl a QDC measurement of the dynamic range of the LED driver is depicted. 
From Monte-Carlo simulations it is known that at the ILC polarimeters signals of up to 
3000 QDC counts are expected in the most extreme detector channels (c.f. Figure H]). 
Thus, Figure [T2b] proofs that the expected dynamic range can completely be covered by 
the LED driver. Furthermore, the dynamic range of the LED driver can be adjusted easily 
with different hlter conhgurations if needed. 

In order to test the stability of the base pulse, the PMT response to light pulses at different 
hxed LED reference voltages has been measured for about one hour each. Figure [13] 
displays two of these measurements at different voltage levels. The data points represent 
the pedestal corrected mean positions of a QDC spectra corresponding to iV = 10® light 
flashes, thus T(x) in the terminology introduced in section [3]2j In Figure [T3al the reference 
voltage has been set from a higher level to Dled = 6.0 V. In contrary, in Figure [TSb] the 
voltage level has been raised from a lower level to Deed = 6.5 V. A dragging behaviour 
in opposite directions can be observed in the hgures although the reference voltage was 
constant. After about 40 minutes the light intensity stabilises. This is important to 
note, since for the proposed nonlinearity measurement the base pulse is altered every few 
minutes. In order to achieve reliable results for the nonlinearity measurement, the base 
pulse should always be scanned in both directions such that the dragging effect can be 
cancelled by averaging. 

The differential pulse stability has been checked on the absolute scale for 12 hours, which 
corresponds roughly to the time of a full nonlinearity measnrement run in the shown 
setup. A drift of at maximum 0.41% has been observed, which is within the tolerable 
range according to Monte-Carlo simulations (see Section 13.2.21) . 
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LED pulse waveform for different LED voltages 




Figure 12: ^ Recorded waveforms of the photomultiplier response for different refer¬ 
ence voltages. The pulses are shorter than 10 ns. (jb]) Dynamic range of the 
LED driver for different hlter conhgurations. The expected dynamic range 
of the ILC polarimeter can be covered. 


Uled=6.5V Uled=7.0V 




Figure 13: Pedestal corrected mean position of the QDC spectra for changing LED 
voltages Dled = 6.5 V and Dled = 7.0 V. About 40 minutes after a voltage 
change a stable mean position is reached. 


During the tests it has been observed that the differential pulse intensity depends on the 
base pulse voltage level. In order to avoid this cross talk for the following nonlinearity 
measurements, two independent LED driver boards have been used. In this modihed 
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setup, only one LED per LED driver board has been used. The other LED on the boards 
has been set to Dled = 0 V. 

Except from the cross-talk problem, the LED driver has proven to fulhl the requirements 
for the usage in the proposed differential nonlinearity measurement. In a future design, 
the two LEDs on one board will be operated in independent circuits. 


4 Measurement Results and Application to the ILC 


In this section, we present a measurement of the nonlinearity for one of the PMTs used 
in the prototype detector, and we discuss the implications for ILC polarimetry. 


4.1 Nonlinearity measurements 

With the setup described in section 13.41 we determined the nonlinearity of one of the 
candidate PMTs to be employed in the Cherenkov detectors of the ILC polarimeters. 
Figure [TTa] shows the result of a differential nonlinearity measurement of the photomulti¬ 
plier R5900-03-M4 operated at 700 V. For each voltage setting of the base LED, responses 
T(a;) to the base pulse alone and T(a; -|- Ax) to base and differential pulse together are 
determined from 10^ samples each. 

As already discussed, a slightly different A T(x) = T(x)— T(x-|-Ax) is obtained depending 
on the scan direction of the base pulse. The black data points show the average of 
both scan directions, which is used to determine the relative nonlinearity in Figure [140 
Since the nonlinearity is dehned with respect to the detector response at the end of the 
dynamic range, the nonlinearity vanishes at this point by construction. It can be observed 
qualitatively that the detector response starts to show a saturating behaviour, which 
is expected for photomultipliers at high light intensities. The measurement has been 
repeated with different hlter conhgurations and exchanged hardware and all measured 
nonlinearities are in agreement with the presented one [25] . 

The measured nonlinearity can now be used for deriving a correction function. In order 
to demonstrate the correction power of this method, an independently measured data set 
has been linearised using the correction function. Thereby, all base-|-differential and base 
pulse measurements have been individually corrected and from this corrected data A T(a;) 
has been calculated again. Figure [T5al shows AT(a;) with respect to the base pulse mean 
position before and after the application of the linearisation algorithm. It is clearly visible 
that after the linearisation the data points for AT(x) are much more consistent with a 
constant than before. Figure ITSbl visualises the originally measured nonlinearity as well 
as the residual nonlinearity after correction, which is more than a factor 5 smaller than 
before the correction. 

This demonstrates that the LED driver, the nonlinearity measurement and the lineari¬ 
sation method are suitable to measure and correct PMT responses at the sub-percent 

level. 
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Figure 14: ^ Measured A T(x) in the default configuration for different scan directions 
of the base pulse, (jb]) Determined relative INL. The absolute scale is not 
accessible with the presented method. 
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Figure 15: ^ Independent data set before (red) and after (green) linearisation. (|b|) 
Visualisation of the measured nonlinearity from Figure [Tib] and the residual 
nonlinearity after the linearisation. Linearity has been improved by about a 
factor of 5. 
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4.2 Application to Polarimetry 

As a final step, we consider the application to the ILC sitnation. For this pnrpose, the 
measnred nonlinearity from hgure Il4bl has been applied within the fast Linear Collider 
Polarimeter Simulation LCPolMC mm- Figure [16] shows the polarisation calculated in 
each detector channel according to equation [T] based on three different estimators for the 
number of Compton electrons N: The true number of Compton electrons (blue, Pc.e.), 
the mean values of the QDC spectra simulated with the measured nonlinearity from 
hgure Il4bl (red, Pmeas) and with the residual nonlinearity after correction as shown in 
hgure [T^ (green, Pcorr)- In addition, the combined polarisation measurement obtained 
from the weighted mean of the per-channel polarisations according to equation [8] is given 
in hgure [T6l 



channel 

Figure 16: Simulation of a polarisation measurement with the ILC upstream polarimeter 
based on the true number of Compton electrons per channel (blue, Pc.e.), the 
mean values of the QDC spectra simulated with nonlinearity (red, Pmeas) and 
after applying the correction (green, Pcorr)- 

In this particular case, the measured nonlinearity causes a deviation of only \6V/V\ = 
0.08%, which would still be within the envisaged uncertainty budget. Comparison with 
hgure [5] shows, however, that the nonlinearity measured for this specihc PMT has a rather 
benign shape, since for an INL of 1% a mean deviation of \5V/V\ = 0.14% is expected, 
with an RMS of ±0.4%. 

Considering the channel-dependency, the measured polarisation values are most obviously 
distorted in the channels with large absolute diherences between the count rates for par¬ 
allel and opposite laser and beam helicities, i.e. at small channel numbers. While at a hrst 
glance it could be argued that channels which measure a polarisation diherent from the 
average could simply be excluded, a closer look reveals that this is not sufficient: For all 
channels with a signihcant statistical weight, and especially those above the zero crossing 
of the asymmetry (channel numbers larger than 10), the uncorrected polarisation (red) 
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is consistently smaller than the valne expected from the nnmber of Compton electrons 
(bine). Fnrthermore, the dependence of the measnred polarisation on the channel nnm¬ 
ber is not a nniqne indicator for the presence or absence of nonlinearities, bnt can also 
be inflnenced by misalignments or tilts between detector and beamline [20] • Thns, the 
differential nonlinearity measnrement is essential for a precise monitoring of the detector 
nonlinearity and its separation from other effects. 

After applying the determined correction to the QDC mean valnes, the obtained channel- 
wise polarisation valnes Pcorr (green) reprodnce the reference valnes from the “trne” nnm- 
bers of Compton electrons nearly perfectly. This lays the fonndation to disentangle effects 
dne to nonlinearities from others cansed e.g. by the above-mentioned misalignments and 
tilts in a real polarimeter system. The relative deviation of the average polarisation from 
the Compton electron-based valne becomes 0.006 %, thns more than a factor ten smaller 
than before, and a factor 16 smaller the effect foreseen in the nncertainty bndget in ta¬ 
ble m This shows that the linearisation based on onr test setnp wonld be fnlly snfficient 
for the ILC polarimeters, and even has the potential to rednce the impact of nonlinearities 
to a negligible level. 

Beyond a carefnl initial characterisation of the photomnltipliers in the lab before the actnal 
installation in the ILC polarimeters, an in-sitn monitoring of the detector calibration is 
clearly reqnired to maintain the achieved precision over many years of ILC operation. 
A possible measnrement scheme at the ILC dnring a rnn is depicted in Fignre HT] It 
makes nse of the pnlsed beam strnctnre at the ILC in which each bnnch train is separated 
in time by 199 ms. The polarimeters are idle in this time snch that this period can 
be ntilised for calibration measnrements. Operating the LEDs at different freqnencies 
/ledi = 2/led2 = 30 kHz would allow for measuring N = 2388 data points ( T(a:); A T(a:)) 
assuming that 80% of the time between two trains can effectively be used. In order to 
reach the demanded precision, N = 3 x 10® data points per base pulse voltage have to 
be acquired (c.f. Section r3.2.2|) . which is possible within less than 5 minutes. Thus, the 
whole dynamic range could be scanned back and forth in 10 steps within less then 1.5 
hours in parallel to a normal ILC physics run. 

In order to make this operation mode possible, small modihcations of the LED driver 
board are necessary: Obviously, it must be synchronised with the ILC clock. Furthermore, 
both LEDs need to be able to be triggered independently, which can be easily realised in 
the current design. This decoupling of the trigger circuit could also solve the observed 
cross talk between the two LEDs. 


5 Conclusions 

In this paper, we presented a calibration system for the Compton polarimeters at the ILC. 
The baseline technology to measure the Compton scattered electrons is a gas Cherenkov 
detector, read-out by single- or multi-anode PMTs. One of the key performance param¬ 
eters is the linearity of such a detector. In order to keep the corresponding contribution 
to the systematic uncertainty of the polarisation measurement below 0.1%, the nonlin- 
earities must be controlled to 0.5 % or better. This is required over the whole dynamic 
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Figure 17: Possible trigger scheme for DNL measurements at the ILC. Between two 
bunch trains 2985 measurements of pedestal (blue), base light pulse (yellow), 
and base+differential (orange) pulse could be performed. This scheme allows 
for operating each LED with a constant frequency between the bunch trains. 
During a traversing bunch train both LEDs must not trigger. 


range of the measurement, which stretches over an order of magnitude in the number of 
Compton scattered electrons per detector channel. The most populated channel has an 
average count rate of 100 — 200 electrons per bunch crossing, depending on the achieved 
laser-electron luminosity. 

At the level of 0.5 %, a calibration via an external absolute reference light source is not 
trivial. The problem of calibrating the linearity of the light source can be circumvented by 
measuring directly the differential nonlinearity with two short (few ns) light pulses: The 
hrst pulse must be tunable to cover the whole dynamic range, but no special precision or 
linearity is required. The second pulse has a hxed size small compared to the maximal 
size of the base pulse. Again, its exact size does not need to be known, but it must be very 
stable, with not more than 0.5 % relative drift during a calibration cycle. The differential 
nonlinearity can then be measured from the response of the detector to the base pulse 
alone and to both pulses together. In case of the application at the ILC polarimeters, 
about 10 scan points for the base pulse with 3 x 10® samples each are sufficient in order 
to achieve the required precision. 

An LED driver providing these pulses has been designed and tested successfully in a setup 
using one of the candidate PMTs to be employed in the gas Cherenkov detectors of the ILC 
Compton polarimeters. In this setup, a correction of the nonlinearities to the per-mille 
level could be achieved. Transferring the result to a simulated polarisation measurement 
shows a negligible impact of the nonlinearity on the polarisation measurement after the 
correction. Finally, we propose a scheme how to collect calibration data in-situ during 
the physics runs of the ILC, exploiting the rather long gaps of 199 ms between the ILC 
bunch trains. A full calibration would then be completed within 1.5 hours, while more 
frequent updates of the calibration can be obtained by a sliding window technique. 

In summary, the performance of the presented calibration system is sufficient to control 
the contribution of detector nonlinearities to the total uncertainty of the polarisation mea¬ 
surement to much better than 0.1 %. This could reduce the total systematic uncertainty 
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of the polarimeter measurements to 8V jV = 0.2 %. In the future, the calibration system 
should be tested with the full prototype gas Cherenkov detector. 
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